Coherent Tunneling of Atoms from Bose-condensed Gases at Finite Temperatures 
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Tunneling of atoms between two trapped Bose-condensed gases at finite temperatures is explored 
using a many-body linear response tunneling formalism similar to that used in superconductors. To 
lowest order, the tunneling currents can be expressed quite generally in terms of the single-particle 
Green's functions of the isolated Bose gases. A coherent first-order tunneling Josephson current 
between two atomic Bose-condensates is found, in addition to coherent and dissipative contributions 
from second-order condensate-noncondensate and noncondensate-noncondensate tunneling. Our 
work is a generalization of Meier and Zwerger, who recently treated tunneling between uniform 
atomic Bose gases. We apply our formalism to the analysis of an out-coupling experiment induced 
by light wave fields, using a simple Bogoliubov-Popov quasiparticle approximation for the trapped 
Bose gas. For tunneling into the vacuum, we recover the results of Japha, Choi, Burnett and 
Band, who recently pointed out the usefulness of studying the spectrum of out-coupled atoms. In 
particular, we show that the small tunneling current of noncondensate atoms from a trapped Bose 
gas has a broad spectrum of energies, with a characteristic structure associated with the Bogoliubov 
Ch ' quasiparticle u 2 and v 2 amplitudes. 
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I. INTRODUCTION 

t> 

The ability to coherently transfer atoms out of a Bose-Einstein condensate (BEC) in a trapped atomic gas is of 
fundamental importance in the creation of coherent matter wave generators, often called "atom lasers" . Theoretical 
interest has developed toward characterizing the out-coupling stage of such an atom laser, with the hope of developing 
■ efficiency criteria for this process. All this work has been done in the zero temperature limit, using either a non- 
interacting gas coupled to a laser mode ]lj or the modified Gross-Pitaevskii dynamics of the condensate order parameter 
0. A complete understanding of the out-coupling process necessarily must take into account the presence of the 
thermal cloud composed of noncondensate atoms. Recently, the role that the noncondensate atoms play in the out- 
coupling process has been discussed by Japha, Choi, Burnett and Band ||, |j. In their work, a "pair-breaking" 
contribution was obtained, in which a correlated pair of atoms take part in the output process, with one output atom 
appearing simultaneously with the creation of an excitation in the trapped Bose gas. 

Out-coupling experiments are only a small class of experiments that can be performed in which a BEC is coupled to 
another gas of atoms (that may or may not be Bose-condensed). The MIT group ^ |) has performed work in which 
a BEC is scattered from an optical Bragg grating. They show that the excitation spectrum of the scattered atoms 
gives a direct experimental measure of the dynamic structure factor of the initial Bose-condensate. Since the detailed 
form of the dynamic structure factor <S(q, uj) is dependent upon the many-body correlations of the Bose-condensate, 
this kind of experiment provides a means of probing the unique correlations induced by a Bose-condensate j?], ^ . 

Despite the large amount of work that has been performed with regard to the coupling of Bose-condensed gases 
with another (possibly Bose-condensed) gas, there still does not exist in the literature a clean formulation of such a 
coupling. Even though Refs. ||, |[ have sparked interest in this issue, the formalism used in these papers does not 
show in a transparent way which correlations present in a Bose-condensed gas can be probed by such experiments. 
Moreover, the case in which both trapped gases are Bose-condensed was not considered. 

In this paper, we set up a general linear response theory approach for treating the coupling of two Bose gases at 
finite temperatures, which gives a formally exact expression for the tunneling current. The general expressions we 
derive — within linear response theory — are independent of the specific approximation used for the many-body 
states of the two uncoupled, trapped Bose gases. In this sense, our results are the analogue of showing how the 
inelastic neutron scattering cross-section is related to the dynamic structure factor @, ^, |io| a nd the general theory of 
tunneling in superconductors [0, |l^| . The work closest to ours is by Meier and Zwerger |13| , who applied a similar 
formalism in discussing coherent Josephson tunneling between two uniform Bose-condensed gases. In Sections II and 
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[II, we give a detailed account of the many-body theory needed to find formal expressions for the tunneling currents 



in terms of single-atom correlation functions. 



'Electronic address: 
t Electronic address: 



luxat@ptiysics.utoronto.ca 



5rifBn@pl1ysics.utoronto.ca 



2 



We apply our theoretical formalism to a class of experiments involving Raman scattering of atoms between trapped 
Bose gases. We calculate the tunneling current expression at finite temperatures using the Bogoliubov-Popov approx- 
imation for the single-particle excitations in the trapped Bose gases O]. We discuss coherent tunneling of atoms 
between two Bose-condensed gases, as well as the out-coupling of atoms from a BEC in a trap to the vacuum (i.e., to 
free particle states). In this latter case, as expected, we show that our more general formalism reproduces the results 
of Refs. H, Q in the appropriate limit. In Section (vi|, we give an explicit calculation of the out-coupling current based 
on treating the noncondensed atoms within a simple local density approximation (LDA), such as used in Refs. |Q |). 

In essence, we show that besides the dominant contribution from tunneling of atoms out of the condensate, one 
finds a much smaller contribution which arises from the atoms tunneling out of the noncondensate. This latter 
contribution, however, is of special interest since it is directly related to the single-particle correlation function of the 
noncondensate atoms. One component has weight proportional to the Bogoliubov amplitude u E N°(E), where N°(E) 
is the Bose-distribution function for an excitation of energy E, and increases with temperature. The other, more 
interesting, component has weight proportional to the Bogoliubov amplitude v E (l + N°(E)). This latter process is 
formally related to the "negative energy pole" of the Bogoliubov spectral density. This so-called "pair-breaking" term 
j| Q is of interest as a unique signature of the unusual correlations induced by a Bose condensate. It corresponds 
to the out-coupling of an atom and the simultaneous creation of a Bogoliubov excitation in the trapped gas. This 
contribution already is present at zero temperature and, in contrast to the tunneling current from the condensate, 
contributes to an output tunneling current of atoms with a wide spectrum of energies. As the temperature increases, 
the weight of the v 2 E term decreases while the u 2 E term starts to build up (we recall that as the condensate density 
goes to zero, u% — ► 1 and v% — >• 0). The analogue of this pair-breaking process was discussed earlier in the context of 
a Bose-condensed gas of excitons [ [l5| . 

In contrast to our work, a proposal has recently been put forward ]l6| to measure these Bogoliubov amplitudes 
using an optical Bragg scattering technique. In this approach, a two-photon Bragg scattering beam is used to excite 
quasiparticles out of the condensate and then a second Bragg beam is used to probe the altered momentum distribution 
of the sample. Recent work at MIT has beautifully verified these predictions [[l7] . Unlike this situation, in the Raman 
scattering experiment the energy absorbed by the atoms imparts a momentum kick in addition to causing a change 
of hyperfine state, allowing the scattered atom to no longer be trapped. 

Our results give another fairly direct method of measuring the Bogoliubov quasiparticle amplitudes. Furthermore, 
since the weight of the out-coupling current is proportional to the Bogoliubov amplitude v 2 , it is significantly amplified 
for larger values of the s-wave scattering length a such as can be achieved by working close to a Feshbach resonance. 
We hope our work further stimulates experimental interest in this kind of out-coupling experiment, even though the 
current from the noncondensate atoms is quite small compared to that from the condensate. 



II. LINEAR RESPONSE THEORY: GENERAL 



Our goal is to calculate the current of atoms tunneling between two atomic gases coupled through a weak external 
field, which can be treated as a time-dependent perturbation to a time-independent many-body Hamiltonian using 
linear response theory (see, for example, Refs. [jll, fL8[), The total Hamiltonian in the Schrodinger representation is 

H(t) = H +e 1 N 1 +s 2 N 2 + V(t), (1) 

where Hq is the Hamiltonian for the many-body systems, e± and £2 are the atomic hyperfine level energies of an atom 
in gas 1 and gas 2 respectively, N±(2) is the atomic number operator for gas 1(2), given by 

N m = J dr$ (2) (r)Vi 1(2) (r). (2) 

The external probe field V(t) is explicitly time-dependent. We assume that the two many-body systems do not 
interact with one another in the absence of V, so that Ho = H10 + -^20- This amounts to assuming that the two gases 
do not overlap, i.e. the coupling term of the kind V'jV'JV^V'i can be neglected. In the absence of V(t), the number of 
atoms in the two systems is fixed by their chemical potentials fix and fi2- 

The tunneling perturbation V{t) couples the two many-body systems together by introducing a mechanism by 
which an atom can tunnel between the two systems. This perturbation has terms which create an atom in one system 
while destroying an atom in the other system, and vice- versa. An example of the kind of experiment described by the 
tunneling perturbation V(t) is the out-coupling phenomenon considered in Refs. |5|, m. In this experiment, an external 
electromagnetic (EM) field causes transitions between two atomic hyperfine levels, one of which is not trapped. We 
label the quantum states corresponding to these hyperfine levels by |1) and |2). Assuming that hyperfine level |1) 
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is lower in energy than hyperfine level |2), energy will be absorbed from (emitted to) the external field facilitating 
the transition |1) — > |2) (|2) — * |1)). In more general experiments, however, the transition can take place between 
two atomic hyperfine levels that are both trapped. An example of such a process is a stimulated Raman transition 
between the magnetic hyperfine levels |1) = \F = l,mp = —1) and |2) = \F = l,mp = 1) (for a detailed discussion, 
see Ref. Il9]). We will return to the specifics of a BEC coupling experiment in Section IV when we calculate the 
tunneling current between two trapped Bose-condensed gases. 

For the moment it suffices to write down the semiclassical form of the tunneling perturbation V arising from the 
electric-dipole and rotating wave approximations || [l^] , 



V(t) 



drl 7e 



'^i(r)$(r)+h.. 



(3) 



Here u> is the effective energy transfer from the laser fields, q is the momentum transfer from the laser fields, and 
7 (taken to be real) is the interaction strength that is associated with the Rabi frequency of the transition between 
atomic states. We have also introduced a factor e vt that models the adiabatic switching on of the interaction at 
t = — oo. In our final results, we take the limit rj — > + . We use units in which H = 1 and ks — 1, unless otherwise 
specified. 

We next outline the formalism for calculating the tunneling current between two Bose gases (resulting from the 
action of V). The state of the combined system evolves in time according to \ip{t)) = U(t) \ip(—oo)), where the 
evolution operator U(t) is unitary and satisfies the equation of motion 



.dtjjt) 



U(t)H(t). 



(4) 



This equation of motion can be formally solved to give U{t) — Texp ^ — i J ^ dt'H (t 1 ) J , where T denotes the 

time-ordering operator. In the Heisenberg picture, operators are written in terms of the corresponding Schrodinger 
picture operators as Ah (i) = U' (t)AU(i), while the state of the system is given by the initial state of the system at 
t = —oo, namely \ip) H = \ip(—oo)). We use the subscript H to denote an operator or state in the Heisenberg picture 
with respect to the full Hamiltonian H . The Heisenberg equation of motion for any operator A is then given as 



dA H (t) 

at 



H H (t),A H {t) 



(•5) 



We define the atomic current from gas I to gas 2 in terms of the thermal equilibrium ensemble average of N%h 

i(t) = -(k 1H (t)). (6) 



We have introduced a negative sign in this expression for the atomic current so that the current flow from gas 1 to 
gas 2 is positive. The ensemble average in (0) is given by 



(N 1H (t)) = Tr{p eg ^ m (t)} 

■[p eq U\t) [#(t),#i] U(t)} 



= iTr 



(7) 



where the initial equilibrium density operator p eq describes the state of the system at t — -co before the coupling V 
is turned on. Since 



H(t), Ni = V(t), Ni , the average in ([?]) takes the form 



(8) 



We now turn to the task of evaluating this ensemble average within linear response theory [ pd| [l8| . 

We first write the average (^J) in an interaction picture, in which V is taken as a perturbation to the many-body 
Hamiltonian H' Q = Hq + e\N\ + £2-^2 of the uncoupled gases. Later we will introduce an interaction picture involving 
K = Hq — HiNi — P2^2- With respect to H' , we define the interaction representation as 



(9) 
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In this interaction picture, we also define the time-dependent density operator as 

Using these definitions, we now can write (||) as 

(^(()) = iTrjp^y^ [v(t),N t ] e"^ 4 } 

Evaluating the commutator V(t),N\ in (|TT|) , we obtain 



N 1H (t))=i / dr 



je^-^Tr { PK {t)4, 1H , (r, t)^ 2H , (r, t)} 



The density operator in the interaction picture defined in (10) satisfies the Liouville equation 

dpm(t) 



dt 



PH>{t),V Wo {t) 



which has a formal solution 

PH> Q (t) = P(-OO) + i I dt' p H , (i'), V H , (f) 

Iterating to first order in V, we obtain 



PH' {t) = Peq + Sp H ' (t). 



Here 



5p H > {t) = i / dt' Pe g ,V H ' (t') 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



is the perturbation to the initial equilibrium density operator p eq to first order in V , the linear response approximation 
to the density matrix. 

Within this first-order approximation for Ph' (t) the ensemble average of an operator A is given by (Ajjit) 



A H , o (t)j +6(A H , o {t)j, where 

(i^(i)) e? = Tz{p eq A H , o {t)) (17) 

s(M(t)) = T r {5p H ,(t)A H ,(t)}. (18) 

The atomic current given by (^]) can thus be written in linear response theory as 

I(t)=I eq (t)+SI(t), (19) 

where I eq is the equilibrium current and SI is the linear response current. Substituting our approximation for PH' {t) 
of (^5|) together with (||) for V into (|l^) , we obtain, after some algebra, explicit expressions for the equilibrium current 

I eq (t) = 21m J d n e^ r -" *> (^urj(r, 4)$^ (*>*)) (20) 

and the linear response current 

£J(t) = 2Re y dr J dr' J dt'e^ t+t '^ 



2 e i[q-(r+r')-a)(t+i )] 



^ (r', tO^Lj (r', 0.^5 ( r , ( r > *) 



(21) 
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The expressions in @ and @ are the basis of the rest of this paper. One sees that there are two different 
contributions to the linear response current in (fjl|), corresponding to phases depending on the times through t — t' 
or t + t' . The former will be seen shortly to correspond to a normal current of single atoms between the two gases, 
while the latter corresponds to an anomalous current of atoms between the two gases. 

We next want to express the ensemble averages in ( p0| ) and (^l|) in terms of thermal Green's functions given by 
standard techniques in finite temperature many-body theory [ fill |18fl . To this end, we first express the time-dependent 



field operators in (20) and (J2l|) in the usual "Heisenberg picture" used in many-body theory, namely the interaction 



picture in which the free Hamiltonian is taken to be the many-body "grand canonical Hamiltonian" 

K a = H' - mNt - p 2 N 2 . (22) 
Using the identity e XNi ipi(r)c~ XNi = c~ x ipi(r), it is easily verified that 

A H ,(v,t)=e-^+^ t iP lKo (r,t), (23) 

where we define ^iK (r,t) = e lKatr ipi{v)e~ lKot . Using (|23|), it is easy to rewrite ( pl| ) in terms of the 4>iK operators. 
We also note that the grand canonical ensemble average of a product of operators of the two gases factorizes 



A 1Ko A 2Ko ) ^(A 1Ko ) <A 2Ko ) , (24) 

/ eq \ I leq \ l 2eq 

where (•••) leg = Tr {pi eg (•••)}. 

Below the transition temperature of a Bose gas, the new phase of matter is described by a symmetry-broken state in 

which the equilibrium density operator does not commute with the number operator: p eq , N ^ (the state breaks 

the system's U(l) symmetry). Thus, the equilibrium ensemble averages are truly over the 77-ensemble introduced by 
Bogoliubov |^(| , which is symbolically defined as 

(Ai Ko ) eq = lim Tr {p v Uq A 1Ko } , (25) 

where 



c 



-/3l(- H 10+-fflai,- i UlJVi) 



Pl q = ^ (26) 

is the symmetry-breaking grand canonical equilibrium density operator. The symmetry-breaking Hamiltonian is 
defined as @ 

H lsb ee J dr (<(r)^(r) + r/ 4 (r)^(r)) . (27) 

Z\ ee Tr-|e~^ 1 ^ 10+ ^ lsi ' _A11 * 1 ) I is the broken symmetry grand canonical ensemble partition function. We will use 

the notation (• • •) to denote equilibrium ensemble averages, with the understanding that they are over a broken 
symmetry state. An equivalent way of stating this is that atomic field operators can be decomposed into two parts: 

^ 1Ko (r, t) = $i (r) + 4> 1Ko (r, t), (28) 

where 

$ (r) = U Ko (r,t)) =^{r)e ie ^ (29) 

is the symmetry-breaking order parameter field and tp is the fluctuation field operator associated with the nonconden- 
sate atoms. In (^), n c0 is the local condensate density while 9 c o is the condensate phase. In our formalism based on 
linear response from thermal equilibrium, the Bose order parameter in ( p8| ) is that appropriate to thermal equilibrium. 
A phase 9 c o(r) would arise if we were dealing, for example, with a condensate with a vortex. We also note that the 
usual time-dependence of the Bose-order parameter (e -IMt ) is absent because we are working in the Kq interaction 
representation. 

We do not need to assume that the inverse temperature is the same for both gases. Indeed, because the two gases 
are assumed to be uncoupled in the absence of the perturbation V, no mechanism exists to force both gases to reach 
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a common thermodynamic equilibrium. Therefore, the equilibrium ensemble averages in (p(i|), ( pl| ) and ( pi] ) may be 
viewed as being over two independent gases. 

The equilibrium tunneling current in (EG) is first order in the perturbation V, whereas the linear response current 



in (21) is second order in V. Using (^3|), (|2C]) can be written as 

I eq (t) = 21m f dne^e-'f"-^-^'' (&x„(r,t)) ($* M) 



Icq 



2eq 



where Ae = ei — e%, A/i = /i 2 — Introducing (|29|) into (30), the equilibrium current is reduced to 



legit) = -27 



J dr^n c io(r)n c2 o(r) sin 



(« - Ae - A/i)< - q r + 6» c20 (r) - clo (r) 



(30) 



(31) 



The magnitude of this equilibrium current is seen to be proportional to the overlap of the two condensate density 
profiles and thus serves as a probe of the condensate surface structure. 

In a non-Bose condensed gas, the ensemble averages (ipiK (r,t)) and (ip2Ko( r ,t)/ m ( p0| ) vanish. As noted 

\ / leg \ / 2eq 

in ( p9| ) and the discussi on p receding it, these ensemble averages do not vanish in a Bose-condensed system. Thus, 
the equilibrium current (|3l| ) is a direct manifestation of Bose broken-symmetry. This expression for the equilibrium 
current is a more general form of the linearized "internal" Josephson current |2l], |2^] for Bose-condensed gases. Our 
expression (|3l]) differs from previous results in several ways. Firstly, there is the phase associated with the momentum 
q transferred by the two laser beams. This will introduce momentum kicks (Bragg scattering) to the atoms in 
the two condensates M. Secondly, we have allowed for the case of a spatially varying condensate phase difference, 
A0(r) = 9 C 2o(r) — 9 c i (r), such as would arise if one of the condensates had a vortex. 

The quantity u> — Ae — A/i in ( |3l| ) is analogous to the vo ltag e in the standard discussion of the Josephson current 
in superconductors Q . We will show later (see ( |59|) and (|60|)) that this quantity also enters into the phase of the 
anomalous component of the linear response current in ( |2l| ) as a second harmonic. We may interpret the atomic energy 
level splitting of the atomic states Ae in addition to the chemical potential difference A/i of the two Bose gases as a 
negative voltage impeding the atomic current flow from gas 1 to gas 2. The energy difference between the laser photons 
acts as a forward biased voltage serving to negate this reverse bias. In the regime where u> — Ae — A/i — > 0, we recover 
what is termed the d.c. Josephson effect in superconducting systems [11 . In this regime, as with superconductors, a 



non-zero difference between the phases of the two condensate order parameters drives a d.c. atomic current. 



III. TUNNELING CURRENTS IN TERMS OF CORRELATION FUNCTIONS 



We next turn to the linear response current in (|2l|), which is second order in the tunneling perturbation V. We 
first want to write it in terms of thermal Green's functions by introducing correlation functions of the atomic field 
operators. As in the previous section, we rewrite our time-dependent interaction picture field operators in terms of a 
Heisenberg picture defined with respect to Kq, rather than H' (see (p3|)). 

We consider the first term in the expression of the linear response current in (^l]). We refer to this term as the 
normal current, for reasons which shall become clear. The normal current is then given by (Ae = e<i— £x, A/x = fi^ — Mi) 



8I N (t) ee 2 7 2 Re J &t'e^ t+t ^ J dr' J dr"e tq < r - r ' 



-i(u 



-A^Xt-i') 



C%°- (r', r, t' - t)C Ko -, (r', r, t! - t) 
-C K % (r, r', * - t')C K ° (r, r', t - t') 



(32) 



We have evaluated the commutators in (^lj) making use of ( p4[ ) to write the ensemble averages as a product of 
two-point correlation functions such as 



C |lV r ^M^^)VWr')) ieg • 



(33) 



In subsequent discussion, we will refer to ipK (r,t) as ip(r,t) for notational simplicity. The normal current given in 
([32"|) is independent of the time t. We make this more explicit by rewriting the integral and noting that factors of e 1 * 4 



7 



factor can be neglected in the limit 77 — > + . Introducing the (Heaviside) step function 9, the normal current may be 
written in the explicitly time-independent form 



/oo r p 

dre-" T e'(r) J dr J d r > e i<i<r-r') e -i("-^-^)-r 



C+tb. ( r '' r ' - r ) C Ml ( r '' r ' _r ) - ( r ' r '' T ) C ^ 2 (r, r', T ) 



(34) 



We observe that, as expected for a current, the second term in ( J34| ) breaks time- reversal symmetry; the current in the 
negative time direction is the negative of that in the positive time direction. 

The correlation functions appearing in (53) may now be written in terms of their Fourier transforms 



/OO J 



(35) 



and hence (|34| ) becomes 

/OO pOO J / /> OO j // PC 

drj ^Lj ^_ e -"^(r) j dr y drV^^V^"-^ 

( r '> r ' "^')%^t (r' 5 r, a/') - CT^t (r, r', u/)C^ 2 (r, r', 



The integral over the time r in (pj]) is trivially given by 

^ T 0r T \ e i{u" —u' -u+As+An+iri)T 



uj" — uj' — uj + Ae + A/i + i?7 



In the physical limit r\ — > + , we may use the identity 



1 ( 1 

lim = P - 

?;^0+ UJ + 17] \LU 



— iit5(uj), 



(36) 



(37) 



(38) 



where P denotes the principal value of an integral. Therefore, the normal current (pq ) may be written in the form 



SI N = SI Na + 81 



Nb, 



(39) 



where we have introduced the two contributions 



/p p OC J / 

dr y dr'e iC i ( r - r ') J — 



and 



C 4l^ ( r '' r ' - w ') C m I r '' r, + ^ - Ae - A/j 



"C^l ( r > r '> w ')^ 2 ( r, r', -a;' - u + As + A/i 



00 da/ •" 



(40) 



2tt 



<5ijv6 = -2 7 2 Im / dr / dr'c^ r - r ')p 1^1^- 



2tt 



( r '< r < ""'^J ( r '' r ' - C «l ( r ' r '' W ')C^ (r, r', -a/') 



w" - w' - w + Ae + A/i 



(41) 



We next show how the correlation functions in (40) and (^lj) may be converted to spectral densities, which in turn 
can be directly calculated from time-ordered single-particle Green's functions jyj. Thus, the problem is reduced to 
the application of well-known techniques from many-body theory to calculate these Green's functions [l8| . We first 
rewrite the correlation functions in terms of the order parameter fields and the single-particle correlation functions of 
the noncondensate field operators. The diagonal correlation functions above become 



C 4t4 (r',r,t) = *S(r')*o(r) + C^(r' ) r,i) 



(42) 
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and 



C u , (r' ; r, t) = $o(r')$o(r) + C^ t (r', r, i). 
Using the detailed balance identity Q 

in addition to the Bose identity 

one may easily write G^ t ^ in terms of a spectral density (see p. 7 of Ref. fl2l|) 

C m (r', r, - W ) = iV ( w )A^(r', r,w). 
Here N°(u>) is the Bose distribution function 

and the single-particle spectral density is 



2 /3w _ 1 



dte 1 ' 



<:<} 



(43) 
(44) 
(45) 

(46) 
(47) 
(48) 



Note that for an inhomogeneous system, this spectral density is in general not real. 

In turn, A^^ is given by the difference of the retarded and advanced Green's functions, G R and G A respectively, 
for the noncondensate quantum field operators of a Bose gas (see p. 150 of Ref. 



A r t r (r', r, u) = i ( G fl (r, r', w + w?) - G A (r, r', u - irj) 



where G R (r, r', uj + 277) is the Fourier transform of 

G>,r',i) =- l ^)([^(r,t),^ t (r') 

and G" 4 (r, r', u> — irj) is the Fourier transform of 

G\r,T f ,t) = iO{-t)(\fc,t),ft{T f ) 



eq 



(49) 



(50) 



(51) 

The point of all this formalism is that standard finite-temperature many-body theory |]l8f most naturally determines 
thermal Green's functions at imaginary Bose Matsubara frequencies iv n = i2nir/f3, n = 0, ±1,±2, . . .. Fortunately, 
one can determine G R ' A (u> ± irj) by a unique analytic continuation G{iv n — > z — > ui ± irj) — G R ' A (ui ± irj) to the real 
frequency axis. This procedure thus gives the spectral density A. 

We may express ( |40| ) and (|4l]) in terms of the condensate order parameters and single-particle noncondensate 
spectral densities. We find 



SI Na = ^Rc J dr J dr'c^-^ 



•■N 



"(r,r',uj — Ae — A/z) 



and 



dw 



/ -nc—nc/ n „/,,/,,/ 



(r, r', u/, w' + cj - Ae - Apt) 



(5/jvb = -2 7 2 Im / dr / dr'c iq ( r - r '^ -!-P 



Alu* t c N nc (r,r',Lu>) 



, 27r a/ - w + Ae + Aix 
2vr cj" - u/ - lj + Ae + Afj, \ ' 



00 j^j' r°° Ai.Ji jnc—nc 



(52) 



(53) 
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In these expressions we have introduced the normal current densities 

\v, r', W ) = $ 10 (r')$io(r)^t^ (r, r', «) - S^r)**,^)^ (r'. r. ). 



N 



(54) 



which corresponds to atomic tunneling of single particles between the condensate and noncondensate components of 
the two gases. In contrast, 



A' 



: (r, r>, u ') = (n?(u) - iV 2 V)) (r', r, u)A^ (r, r', c/), 



(55) 



is associated with atomic tunneling of single atoms between the noncondensates of the two gases. 

We observe that in both <5//v a and SI^ b , atoms make only two types of transitions between the states of the 
two gases. The first transition process involves tunneling of an atom between the condensate of one gas and the 



noncondensate of the other. It is represented in 5In ci and <51jvb by terms of the form ^*, 2 ^ 1 ^ 2 )Aj l 



^2(1)^2(1)' 



This 



contributes at all temperatures including T = 0. The second transition process involves the tunneling of an atom 
between the noncondensates of the two gases. It appears in 5If^ a and Sl^b as a term of the form ^A^t This 
kind of transition can occur at zero temperature only when one of the gases is Bose-condensed. 

Finally we consider the "anomalous current" given by the second term in ( pl| ) in terms of the condensate order 
parameters and noncondensate spectral density functions. We follow the identical procedure used to find ( [52] ) and 
([53|) for the contributions 51 n a and Sl^b to the normal current. We obtain finally 

5I A {t) = 2 7 2 Re J dr J dr'c lq ( r+r ') J dt' e -*'(«>- As- A„+»7) (*+*') 

x C^{v',v,t' -t)C mi {v'^t' -t) 

-%(^',i-('%»(r,r',(-f')]. (56) 

Clearly this current may be written entirely in terms of the so-called anomalous correlation functions involving creation 
or destruction of pairs of atoms, hence the name anomalous current. Using t + t' = 2t + t' — t = 2t + r, we first rewrite 

SI A {t) = 2 7 2 Reje^ 2( "- A£ - /Vl)t y dr J drV q ( r+r,) J dr6»(r)e l( "- A£ ~ A ' J+l ' 7)T 

x C Ml (r', r, -r)C m (r', r, -r) - C Mi (r, r', r)C^t (r, r', r)] | . (57) 

This shows that the anomalous current is explicitly time- dependent, in contrast to the normal quasiparticle current 
given by ([36]) . It clearly involves a second harmonic of the first-order Josephson current in (|Tj) . 

Following the same procedure we used in writing the normal current in terms of the condensate order parameters 
and noncondensate spectral density functions of the of the two gases, the anomalous current may be rewritten as 



6I A {t) = 6I Aa (t) + SI Ab (t), 

where we have introduced the quantities 

5I Aa (t) = 2 7 2 J dr J dr'j Aa (r, r', ui — Ae — A/j,) 



(58) 



2(w 



A/j,)t - q ■ (r + r') 



(59) 



and 



5I Ab (t) = 2~f 2 J dr J dr'J Ab (T,v' ,lu- Ae - A M ) 

2(w - Ae - Afi)t - q • (r + r') 



x sin 



(60) 
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The quantities J Aa and J Ab in the above expressions for 51 Aa and SI Ab are given by 

1 f°° duj' 

J Aa (r,r',uj- Ae- A/x) = - / — Reu(r,r»'-u + A £ + A/i) 

^ J-oo 27T 

00 du/ f 00 duj" lmi A (r,r',uj',u") 



2tt 2?r uj" - uj' + uj - Ae - A/i 



(61) 



(62) 



and 

1 f°° dec/ 

J A b{r,r',uj) = - / — Imi A (r, r', w', w' - w + Ae + Ap) 

2 J-oo 27r 

+ J_ oc 2vr J_ oc 2ir uj" - uj' + to - Ae - A/x ' 
where we have introduced the total anomalous current density iA, which is written as 
u(r,r',w,o/) = ^w(r)^io(r')A^(r,r',oj')8{oj) 
-$l o (r)$| o (r')%^(r',r, W )5(o;0 

+ (^V) - ^V)) ^ ( r '> r > <») A Ml ( r > r '> <"')• (63) 

Here A^(r',r,uj) is the anomalous noncondensate single-particle spectral density function, analogous to the noncon- 
densate single-particle spectral density function A^^ defined in ([i"8|). 

The first and second terms in (63) correspond to a tunneling current between the condensate of one gas and the 
noncondensate of the other gas. The third term in (^) represents a tunneling current between the noncondensates 
of the two gases. 

One can rewrite J Aa and J Ab such that the contributions from the condensate-noncondcnsate and noncondensate- 
noncondensate tunneling processes are made more explicit. Thus, we have (see also Ref. ]l3|]) 

1 t,„,c-»c_ _/ A . InT dw' Imzr ic (r,r>') 



J Aa (r,r>-Ae-A//) = i- Re^-" c (r,r',^ - Ae - A/x) + 



27r w' - w + Ae + A/x 



i r 00 dw' 



+i / ^Re l ^" c (r,r',^,c'- W + Ae + A M ) 

2 J-oo 27T 

duj' [°° duj" Im^ c -" c (r,r>'V) 



and 



2?r 2tt cj" - w' + uj - Ae - A/i 



(64) 



1 1 /' c 

J A i(r,r', w -Ae- A/x) = — Imx^" c (r, r', u - Ae - A/x) - — P / 



27r w' - uj + Ae + A/x 
+ - / — Im z^ c "" c (r, r', uj' , uj' - uj + Ae + A/x) 
cV Z" 00 dw" R ei X"™ c (r,r',^',^") 



(65) 

2tt 2?r w" - uj 1 + uj - Ae - A/x' v ; 

where we have introduced the condensate-noncondensate current density (the analogue of ([54])) 

x^" c (r, r', W ) = $ 10 (r)$ 10 (r')^t^ (r, r', - w ) - ^(r)^^)^^ (r', r, w ), (66) 
and the noncondensate-noncondensate current density (the analogue of ml 



x^- c (r, r',uj, uj') = (n°(u) - xV 2 V)) A^ (r', r, uj)A^ (r, r',uj'). (67) 



We note that SlAb in (|60j) is the current most directly associated with the coherent tunneling responsible for the 
Josephson effect in superconductors. This current is in-phase with the Josephson current in (|3l]) but involves a second 
harmonic with a frequency 2(uj — Ae — A/x). It clarifies the role of the noncondensate in the BEC Josephson effect. 
On the other hand, the SIah current in (B9) is directly related to a dissipative interference current, just as one finds 



for the Josephson effect in superconductors 112, 13 
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IV. COUPLING BETWEEN TWO TRAPPED GASES 



In this section, we use the results of Section [II to calculate the second order tunneling current involved in a Raman 
scattering experiment. The spectral density functions are calculated in the finite temperature Bogoliubov-Popov 
approximation pi] . The Josephson tunneling current has already been discussed and is given by ( |3l| ) . 

A two-photon Bragg scattering experiment involves two laser beams at different frequencies and with different 
propagation vectors [pi. However, we are interested in a different type of two-photon process in which an atom 
changes its hyperfine state. This is a Raman process and is depicted in Fig. [j]. We now discuss how the effective 
frequency u> and wavevector q in (J3J) are defined in such a Raman scattering experiment fl9fl . Assuming that each 
laser beam has equal intensity /, the combined intensity of the two beams at space-time point (r, t) is given by 



I tot {r,t) = 2/ + 2/cos(q-r-cj<), 
where uj is the frequency difference between the two laser beams 

The wavevector difference between the two laser beams is 

q = k a k b , 

where k, is the wavevector of laser beam j. We may take k\ = k\ 
relation (see, for example, p. 103 of Ref. pj) holds and q is given by 



(68) 



(69) 



q 2 = 4fc 2 sin 



(70) 

k 2 such that the standard Bragg scattering 



(71) 



where 9 is the angle between the two counter-propagating beams. Finally the coupling parameter 7 in (^|) is related 
to the single-photon Rabi frequency ^luabi (assumed to be the same for both hyperfine level transitions) and the laser 
beam detuning A from the intermediate level |3) (see Fig. P according to |19j 



|^,Rabi| (2) 

7 oc ex ir ' 



Rabi' 



(72) 



(2) 

where ^^bi is the so-called two-photon "Rabi frequency" . 

The many-body Hamiltonian (in the absence of a tunneling perturbation) is Hq = Hiq + H20 , where 



flio = / dr\4(r) 



2m 



+ U ext {r) 



Mr) 



9 ^\{v)^\{v)Mv)Mv) 



(73) 



In this expression for H 10 , U ext is the trapping potential and g\ = Anai/m is the s-wave scattering pseudo-potential 
with ai being the s-wave scattering length for atoms in gas 1. The expression for H20 has an analogous form. We recall 
that the spectral densities in Section III arc defined in terms of the grand canonical Hamiltonian Km = Hio — (J,{Ni. 
Thus the excitation energies of gas i are measured relative to the chemical potential 



A. The Bogoliubov-Popov approximation 



In the standard static Bogoliubov-Popov approximation for excitations in a trapped Bose gas |T 
noncondensate field operator tp can be written as 

^(r, t) = £ (uMe-^aj - v*{v)e^a]) , (74) 
3 

where Uj(r) and Vj(r) arc the Bogoliubov quasiparticlc amplitudes, Ej is the energy of Bogoliubov state j and ctj 
(ctj) annihilate (create) a quasiparticle in state j, respectively. These annihilation and creation operators are defined 
to satisfy the boson commutation relations 

[a m ,a n ] = [d4„,o4] =0; [a m ,o4] = S mn . (75) 
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The Bogoliubov quasiparticle amplitudes must satisfy the finite temperature coupled Bogoliubov-Popov equations of 



motion 14 



Cuj(v) -gn c (r)vj(r) = EjU^r) 
Cvj(r) - gn c (T) Uj (r) = -E 3 v 3 (r) 

in order for Kq to be reduced to the diagonal form 

K B °a = cons t + ^ Ej&t&j. 



(76) 



(77) 



In (J7g), we have introduced the differential operator 



V 2 



£ = --z \- U ext (r) - u + 2gn(r), 

2m 



(78) 



where n(r) is the total static density. 

Using the commutation relations (|75|), one finds that the single-particle spectral density function can be expressed 
in terms of Bogoliubov quasiparticles as 



A m (r',r,t-t') = J2 

3 

which has the Fourier transform 



j 



<(r')u i (r)e- i - B ^ t -*') - u J (r>?(r)e iB '< 1 -*'> 



u*{v')u 3 (v)5{uo - Ei) - Vj (i>)vXr)6(w + Ej) 



(79) 



(80) 



We note the appearance in (|80|) of negative energy poles (to — —Ej) in addition to the usual positive energy poles 
(u> = Ej). These are a result of unique processes induced by the Bose condensate. The two poles follow immediately 
from <pm) , which shows how the destruction of an atom involves a linear superposition of the destruction (amplitude 
u) and creation (amplitude v) of Bogoliubov excitations. 

Similarly, the off-diagonal or anomalous single-particle spectral density function of relevance to the anomalous 
current densities nc in d66| ) and i n £~ nc in ( |67[ ) is given by 



^(r',r )W ) = -27r^; 



h {v)v*{v')8{uj - E 3 ) - Uj (v')v*(r)S(w + E 3 ) 



(81) 



The spectral density functions can be shown to satisfy A^^ (r', r, w) — At-(r, r',u) and A^^(r',r,ui) 



B. The normal and anomalous currents in the Bogoliubov-Popov approximation 

Using the results in ( |so| ) and ( [Bib, w e can now write the normal and anomalous linear response currents in terms of 
Bogoliubov excitations. In SectionUIlL we found that two distinct physical processes contribute to the normal current. 
The current density of atoms tunneling between the condensates and noncondensates of the two gases is given by 
©.or 



$* 10 (v)<t> w (r')u* 2j (r)u 2j (T')5(w - E 2j ) 



-$J (r')$ 10 (r)^(r')^(r)5(a; + E 2j ) 
-*5 (r / )* ao (r)«^.(r')u ai (r)*(w + E l3 ) 

+^ (r)$ 20 (r')v lj (r)v* lj (r')6(u - E l3 ) 



(82) 
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This type of process involves both condensate atoms and quasiparticles. The current density associated with the 
tunneling of atoms between the noncondensates of the two gases has a more complicated form given by (|55|), or 



N 



3 (r,r>, W ') - (27r) 2 (V°H-Ar 2 V) 



E 



u; i (r , )«i i (r)u^(r)u^(r')«(w - - E 2j ) 



-uUv^ui^v^v^ir'Wu - + E 2j ) 

-^ li (r / )^ i (r)u^(r)M 2j -(r / )<5(a; + £«)<V - By) 

+«i i (i>J i (r)v2 j (r)v^(r')«(w + £ u )5(o;' + E 2j ) 



As in the case of the normal current, two different physical processes contribute to the anomalous current in (|5^ 
The first process involves pairs of condensate atoms tunneling to noncondensate states of the other Bose gas. This 
anomalous current density is given by (pq), namely 



(83) 



$io(r / )*io(r)«!w(r)«3 J .(r / )«(w + E 2j ) 



-$l (r')$* w (T)u* 2j (T')v 2j (r)6(uj - E 2j ] 



(84) 



The expression for the anomalous current density between condensate 1 and noncondensate 2 has a similar form. The 
final process involves quasiparticles of both gases and is given by (p^) , namely 



i" c -" c (r,r>,a/) = (2tt) 2 JV°( W ) - JV 2 °(u/) 



E 



un(r / )»I i (r)u^(r')w 2i (r)*(w + E xi )5{J + E 2j ) 



-u u (r')vUr)u 2j (r)v* 2] (r')d(uj + E U )5{J - E 2j ) 
-^ i (r) Uli (r')u* J .(r , )^(r)5(w - E u )5{u' + E 2j ) 

+u* ll {v)vu{v')u 2] {v)v* 23 {v')5{uj - E u )5{lo' - E 2j ) 



(85) 



As with the Josephson current given in (|T]), one sees from (|82|)-(p5|) that the magnitude of the normal and anomalous 
quasiparticle currents is determined by the spatial overlap of the Bogoliubov wavefunctions described by u%{r) and 
Vi(r) for the two systems. 

When one uses (|82|)-(|85|) in the expressions for the normal (see ( p2| ) and (p3(j) and anomalous (see (|5^) and (|6£ 
currents, our results are essentially equivalent to expressions discussed in Rcf. |l3| . The main differences are that 
Meier and Zwerger approximated the two Bose gases as being uniform (or homogeneous) and used the single-particle 
spectral densities at T = 0. 



V. OUT-COUPLING TO A NON-BOSE-CONDENSED GAS 

In this section, we consider the special case of atoms out-coupling to the vacuum as discussed in Refs. j| [|. More 
precisely, we first deal with the case when gas 2 is a uniform non-interacting gas. The corresponding spectral density 
function then has the well-known form jl8[ 



A^ 2 ( k > w ) = 2nS(u;-e k2 ), 



(86) 
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where 

e k2 = k 2 /2m - fi 2 (87) 

is the the excitation energy of an atom in gas 2. Because the gas is assumed to be in the classical limit, we have 
iV 2 °(e k2 ) « 1. 

The only terms in (^2|) which are finite reduce to 

= $io(r')<Mr)A^ 2 (r-r', W ) 

/dk / 
__ e -ik.(r-r )$* ( r ' )(i , lo(r)( 5( £ ~ k2 _ W ) (88 ) 

and in (|§) 

J --r(r,r',^c') = N^)A m y,r^)A m2 {v-T',J) 

j 

+ [1 + iV"^)]^^)^ ^)^ (w + £y) W 2 - w'). (89) 



We have used the Bose identity 7V°(-cj) = -[1 + N°(u)] in writing the last term in (||). The condition that iVf < 1 
results in negligible tunneling from gas 2 back to gas 1 . 

The current density i^ ^ in ( j§8| ) corresponds to the destruction of a condensate atom in gas 1 and the creation of 
an excitation in gas 2 with energy £k2- The first term in the density ^out of (|39|) represents what might be called the 



Nout 

"quantum evaporation" of an atom, destroying a Bogoliubov quasiparticle with energy Eij = k 2 /2m— (uj — Ae+fii) and 



creating an atom in gas 2 with k /2m. The second term in (89) corresponds to creating an atom in gas 2 with energy 
k 2 /2m and simultaneously creating a quasiparticle excitation in gas 1 with energy Eij = — [k 2 /2m — (u — Ae + fix)]. 
This second out-coupling channel in ( |89[ ) is a direct consequence of the unusual condensate-induced correlations in 
a Bose gas. It corresponds to the "pair-breaking" process discussed by the Oxford group ||. The analogue of this 
process appears in the predicted lineshape of the recombination photons emitted in the decay of an exciton in a Bose- 
condensed gas of optically-excited excitons ]lq|. All three out-coupling contributions will be calculated in Section 



VI 



Using (p8p and (|89|), the out-coupling current 5 1 out is 

Slout = 7 2 /dr / dr' / -^e^M'-'' 



(2tt) 2 

j$ 10 (r')$io(r)<5 (k 2 /2m - (u - As + 
+2tt^ A^ijK^rOwi^r^ (fc 2 /2m - (w - Ae + Mi) - Eij) 

3 

+ [1 + iV^SuOlfi^rO^W^ (k 2 /2m -{u;-Ae + /xi) + By) 



(90) 



We refer to the second term as the "u 2 contribution" and the third term as the "v 2 contribution" . It is useful to 



describe the main features of the out-coupling current spectrum arising from the three terms in (90), as a function 
of the energy k 2 /2m of the out-coupled atoms, for fixed u> and q. The energy conserving delta functions in (|9p| ) 
determine the energies of the out-coupling current peaks for the different processes. There is a delta function peak 
at k 2 /2m + £2 = ^ + £1+^1 from atoms coming directly from the condensate, where \x\ is the condensate energy 
of a trapped atom in the hyperfine state |1) and e% is the hyperfine state energy. Relative to this condensate pole, 
the peak associated with the quantum evaporation term is shifted by a positive amount Eij, while the peak of the 
"pair-breaking" channel is shifted by a negative amount —Eij. These shifts follow from the fact that the quantum 
evaporation channel arises from the positive energy pole (related to u 2 ) and the pair-breaking channel corresponds to 
the negative energy pole (related to v 2 ) in the single-particle spectral density. 

We next consider the behavior of the three out-coupling current contributions in ( |90| ) as a function of temperature. 
For T\ = 0, since N\(E\j) — > 0, (i.e. no thermally excited quasiparticles) the second channel (u 2 ) has zero weight. 
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Thus, at Ti = 0, only the first and third processes occur. When gas 1 is at a finite temperature (Tj ^ 0), the second 
channel (u 2 ) begins to contribute to the out-coupling current, taking weight away from the pair-breaking channel 
(v 2 ). All these features are shown explicitly in the numerical results given in Section VI (see Figs. 

The out-coupling current derived in Ref. Q is given as the rate of atoms excited out of gas 1 with energy £k2 and 
may be written as (where we use a notation more consistent with our own) 

^ = 7 2 /dr| dr'e^(- r ')^ 2 (r)^ k2 (r') 

x |$ 10 (r')$io(r)<5 (k 2 /2m - (w - Ae + M i)) 
+2ir^2 JV?(J5y)«J j (r')ttij(r)5 (k 2 /2m -{uj-Ae + m) - E Xj ) 

3 

+ [l + N°(E l3 )]v l3 {r')v* l3 (r)6(k 2 /2m-(LU-Ae + Li 1 ) + E l3 ) j, (91) 

where (/?k2(r) is the wavefunction for the out-coupled excitation with energy £k2- We can make a direct comparison 
between the integrands of ( |90| ) and (|9l]). Indeed taking y>2k(r) = e jk r (a free atom with momentum k) reduces ( pl| ) 
to the integrand of (p0|). Thus our result for the out-coupling current reduces to that derived in Refs. || ||, where 
the out-coupled atoms were in free-particle states. 



VI. THE LOCAL DENSITY APPROXIMATION (LDA) 

In this section, we evaluate the out-coupling current in ( |90| ) and estimate the magnitude and temperature- 
dependence of the quantum evaporation and pair-breaking effects discussed in Section In previous work H g, these 
noncondensate processes were found to be a few percent, relative to the dominant contribution of atoms tunneling 
out of the high density condensate. It has also been found [ p6| that in current Bose gas experiments, the contribution 
of the noncondensate atoms to the dynamic structure factor S(q, u>) at T = 0.9Tbec is about 5 % of the contribution 
from the condensate fluctuations (the latter has been studied in Refs. || EL || at T — 0). This small amplitude is, 
of course, a direct reflection of the strongly peaked condensate at the centre of the trap compared to the broad, low 
density profile of the noncondensate atoms. However, while small, the tunneling current from the noncondensate 
atoms gives direct information about the Bose-induced correlations and thus is a worthwhile goal for future Bragg 
scattering and Raman experiments. 

The condensate current in (BOft can be rewritten in the form 



5I c ou M,5)^l 2 J ^|$ 10 (q-k)| 2 5K), (92) 

where $io(q — k) is the Fourier transform of the Bose order parameter 

$io(q - k) = J drc^- k ) r $ 10 (r) (93) 

and we have introduced the abbreviation = k 2 /2m — (u> — Ae) — ^l\tf- It is also useful to introduce a "detuned 
frequency" 

8 = uj - Ae. (94) 

The noncondensate contribution in ( po|) is more conveniently written in terms of the trapped gas Bogoliubov spectral 
density in ( |80| ) , 

61^,5) = 7 2 | dr J dr'| ^e^-W^N^A^r',^^). (95) 

In the rest of this section, for notational simplicity, we shall drop the index "1" on the quantities associated with 
the trapped gas. Equations (p2\) and (95) summarize our main results for the tunneling currents into the vacuum, 



showing very clearly which properties of the trapped Bose gas they depend on. 
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A. Condensate contribution 

In the standard Thomas-Fermi approximation, the condensate density profile is given by pa 



where the chemical potential is 



and 



n c0 (r) = < V MTF " Uextir) ) /9 > T - RtF , (96) 
0, r > R TF 



fi TF = ^muj^R TF , (97) 



x V5 

a ■ 



Rtf = a HO 15JV C (T) . (98) 

V a-HO J 

The Thomas-Fermi radius Rtf defines the size of the condensate, with N C (T) the total number of particles in the 
Bose-condensate at temperature T and ano — y/h/mtjjQ is the oscillator length. The T = Thomas-Fermi results 
have been extended to finite temperature by adjusting iV c according to the noninteracting trapped gas result (25) 

Nc(T) -l-(J-)\ (99) 



N \Tbec, 

where ksTBEC = 0.94ftwo^ 1 ^ 3 is the noninteracting trapped gas Bose-condensation temperature. Using ( |96|) for the 
order parameter $o(r) = y/n c o(r), the Fourier transform in (^3|) can be calculated analytically 

-sin(fcr) 



$o(k) = 4tt / drrVra c0 (r)- 



(i 



2?r2jR ^j 2 (fci? TF ), (100) 



where Ji is the ordinary Bessel function of order 2. This result is the same as previously obtained J8|, p7| , taking into 
account that we use harmonic oscillator units (length is measured in units of clho = yh/rnujQ and energy is measured 
in terms of fajo). 



Writing the integral in (92) in terms of spherical coordinates, we obtain using (jlOCj) 

7 2 ( 5 n c0 (r = 0)) 7/2 



M5ut(q,<*) 



2-kuq g 

x- / dx^-6(6 + flTF), (101) 
q J x _ x 4 



where x± = y / 2gn c o(r = 0) [^5' z ' = + A*tf) ) • The result (101) may be easily generalized to asymmetric harmonic 
traps, having trapping potential 

U ext (r) = im^ 2 + W y+^ 2 j, (102) 
by introducing the coordinate transformation 

x = — x' , y — — y' , z — — z', (103) 

LU X LU y LU Z 

where u> = (uj x LUyUJ z ) 1 / 3 ' . The same procedure as used above follows through with the same final expressions, apart 
from ujq being replaced by Q. In addition, any wavevector k in our tunneling expressions such as ( |92"| ) needs to be 

interpreted as k — > (jj-k x , ^k y , ^-k^j in an anisotropic trap ||. 
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The total out-coupling current from the condensate given by (j92|), as a function of the detuning parameter 5 in 
(|94|), is shown by the solid lines in Figs. || and ||. In all our numerical calculations, for the two-photon Rabi coupling 
we have used 7 = 0.2u>o. Note the logarithmic scale for the out-coupled atom current. The current magnitude only 
decreases as a function of increasing temperature by about 20%. From Fig. ^ it can be seen that the maximum 
condensate tunneling current amplitude is sensitive to the magnitude of the momentum transfer q. The decrease in 
the maximum tunneling current with increasing values of q is a result of the fact that the momentum distribution of 
atoms in the condensate is peaked at zero momentum, decreasing for increasing momentum. For an asymmetric trap, 
this means that the directions of strongest confinement (largest trap frequency) will out-couple most efficiently. We 
observe also that the peak of the condensate tunneling current shifts toward positive S as q increases. 



B. Noncondensate contribution 



The out-coupling current ( p5| ) from the noncondensate atoms will be estimated using a simple local density approx- 
imation (LDA). The LDA has also been recently used to calculate 5(q, ui) in trapped gases |7|, [| . We first introduce 
a new coordinate system defined in terms of the relative coordinate 



r = r — r 



and the center-of-mass coordinate 



R 



r + r' 



(104) 



(105) 



Writing the single-particle spectral density function Ax^x(r, r', ui) in ( p5| ) in terms of these new coordinates, we note 
that for large condensates and quasiparticle energies much larger than Hluq, the spectral density will vary slowly as a 
function of R. It is thus plausible to approximate A^^ by its homogeneous form, but adjusted so that at each point 
R the appropriate condensate density rt c o(R) is used. Thus, in the LDA, the single-particle spectral density function 
in (95) is approximated by 



AM(r,r» = j4 w (f,R,w) ~ A$£(f , u; n c0 (R)) 



(106) 



where A 1 -®-^ denotes the uniform gas expression for the spectral density function for a condensate density n C Q. 



In the Bogoliubov-Popov approximation, the LDA form ( 106 ) for the spectral density in energy-momentum space 
is M 



A^ A (p,Lu;n c0 (R)) = 2vr 



(p;n c0 (R))<5 u;-£(p;n c0 (R)) 



v 2 (p; n c0 (R))S uj + E(p; n c0 (R)) 



Thus, in the LDA the quasiparticle energies become 



£(p;n c0 (R)) 



e| + 2Me//(R)£ P , R< Rtf 

£ P _ Me//(R)i R > Rtf 

while the Bogoliubov quasiparticle mode amplitudes are given by 

f gp-Pe//(R) j_ 1 



(107) 



(108) 



U^(p;n c0 (R)) = <^ 2E(p;n c0 (R)) ' 2 



1. 



^(p;n c0 (R)) = { 2£ 



R 5: Rtf 
R > Rtf 



2£(p;r lc0 (R)) 2> U - K TF 



R > Rtf 



We have introduced the free atom energy e p = p 2 /2m and 

jti e //(R) = fi T F - U ext {R) - 2gn c0 (R), 



(109) 
(110) 

(111) 
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plays the role of an "effective" local chemical potential, modifying the chemical potential \itf as a result of the 
inhomogeneity of the trap. 

Within the LDA, the noncondensate contribution ( |95| ) to the out-coupling current is given by 

where as before we have introduced the variable = £k — <5 — /itf for notational simplicity. Using the identity 

*(/(aO)=E jTT^C* -*0. (H3) 



where f(xi) — 0, we may evaluate the integral over R in ( |112| ) assuming a spherically symmetric trap. After some 
algebra, one obtains the following expression for the noncondensate contribution to the out-coupling current 

2tt 7 2 f dk - LDA 



(57- (q, S) = ^J N^)A^(d - k, uu), (114) 



where we have introduced the LDA spatially averaged spectral density function in energy-momentum space 
Aj^(k, W ) = |dR^(k,u;;n o(R)) 



el+ui 2 



2ttW 2 Utf - /i(k, w) te 2 + M (k, w ) mtf - M(k, w) 



24 £k 



-2nJ2^fi TF - v(k,Lu)J ^- ^^ + M)fi{-^)H[ /,ik.u.-) (//,-,, - //(k.^i 



+47^/2(0; - £ k + n TF )0{ui - e k ), (115) 

where 

M (k,c)^^=-^ (116) 



and harmonic oscillator units are used. Note that the noncondensate out-coupling current in (114) may be generalized 



to an asymmetric trap in the same manner as discussed earlier for the condensate out-coupling current (see remarks 



below (103)) 



The formula in (114) has a simple physica l int erpretation. The integrand gives the number of emitted atoms with 



momentum k per second. The first term in ( 115 ) is associated with the pos itive energy pole in (p0[), with amplitude 



and only occurs for positive frequency uj > 0. The second term in (115) arises from the negative energy pole in 
(pOf), with amplitude v 2 , and only occurs for negative frequency uj < 0. These first two terms are associated with 
quasiparticles in the region of the (Thomas-Fermi) condensate, i.e. for R < Rtf- The third term, in contrast, is 
associated with quasiparticles lying outside of the condensate region. These quasiparticles comprise a Hartree-Fock 



gas in our simple model calculation. The dominant weight in (115) comes from small e k and u>. The quasiparticles 
with u 2 and v 2 weights in the region where there is a condensate are the dominant contribution. For u> > £k(> 0), 
one also finds a contribution from the normal Hartree-Fock gas of atoms existing outside of the condensate. As the 
temperature increases toward Tbec-, this contribution from the normal Hartree-Fock gas becomes more significant 
relative to the Bogoliubov quasiparticles. Of c ours e, for T > Tbec, it is the sole contribution. 

We have evaluated the integral appearing in ( h_14| ) numerically. Some results, as a function of the detuning parameter 
S for various temperatures, are shown in Figs.^j and [| Note the logarithmic scale for th e ou t-coup led a tom current. 



As a function of detuning, the out-coupling current separates into two regions. From (114) and (115), we observe 
that for 5 < —\ltf (or u> + e\ + {Itf < £2), we have v k > for all values of the wavevector k. In this region, we 
find a contribution from the u 2 Bogoliubov modes, in addition to a smaller contribution from the atoms coming from 
outside the condensate region, at finite temperatures. This contribution to the out-coupling current Sl^f represents 
the "quantum evaporation" process discussed at the end of Section |v|, being associated with positive energy poles 
of the single-particle noncondensate Green's function. It is only present at finite temperatures. In Fig. |^, we have 
plotted the total out-coupling current as a function of the detuning d at an intermediate temperature of T = 0.2Tbec- 
However, for S > — (J-tf (or u> + e% + \itf > £2), can be less than zero for a restricted set of wavevectors k. In 
this case, there are two contributions to the out-coupling current SIpf. One arises from the v 2 Bogoliubov modes and 
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hence represents the so-called pair-breaking process associated with the negative energy poles of the single-particle 
noncondensate spectral density function (i.e., a tunneling atom creates a Bogoliubov excitation in the trapped gas). 
This is the only process which contributes to the noncondensate out-coupling current at zero temperature. At finite 
temperature, a second process arises which is associated with quantum evaporation from high momentum states. 
This involves the positive energy poles of the single-particle noncondensate spectral density function. This second 
contribution decays more quickly as the detuning is increased, compared with the v 2 contribution. We note that these 
two regions will be Doppler shifted toward positive <5 by the nonzero momentum kick q from the laser fields. The 
overall behavior of our LDA results in Figs. §-[|is similar to that reported in Refs. j| Q], which were based on solving 
a one-dimensional Gross-Pitaevskii equation for an isotropic trap. 

We have also evaluated the noncondensate contribution to the out-coupling current for q = 10/a#o, with the 
results shown in Figs. ^ and We observe that the point at which there is no longer any contribution from the v 2 
Bogoliubov modes shifts increasingly toward positive 5 as q is increased. However, unlike for the case of q = 2/ano, 
the contribution at finite temperatures from the u 2 Bogoliubov modes dominates that from the v 2 Bogoliubov modes 
for the range of 5 considered. This results from the fact that the higher momentum quasiparticles have greater spectral 
weight. We have plotted the total out-coupling current at a temperature of T = 0.2Tbec in Fig. ||. 

Because of the strength of the condensate contribution, some comment is required about the observability of the 
noncondensate contribution in the total out-coupling current. Since the contribution from the condensate is smallest 
at high q, we compare the condensate and noncondensate contributions at a momentum kick of q = IQ/auo (see 
Fig. |f|). Since the out-coupling current as a function of 5 is much broader for the noncondensate contribution than 
the condensate contribution, there are regions in which the strong contribution from the condensate will not dwarf 
that from the noncondensate. At T = Q.QTbeCi the maximum contribution from the quantum evaporation process is 
about 10% of the maximum contribution from the condensate (see Fig. ||). The small but interesting v 2 Bogoliubov 
contribution is perhaps more visible as a "distinct contribution" in the limit of very low temperatures (see the T = 
results in Figs. || and ^). 

The unique features associated with the out-coupled atoms coming from the noncondensate can also be enhanced 
by working with larger values of the s-wave scattering length a. In particular, one can "tune" the scattering length 
a significantly by the use of Feshbach resonances p9| , B^ |. The amplification of the out-coupling current of atoms 
associated with the negative energy poles of the single- particle noncondensate Green's function could be looked for in 
trapped gases near a Feshbach resonance. Out-coupling experiments may thus provide important information about 
the Bose-condensed gas excitation spectrum if we can increase the s-wave scattering length substantially. In this case, 
of course, we expect observable changes from the simple Bogoliubov-Popov approximation we introduced in Section IV 
for the single-particle spectral density. Such many-body effects can thus be probed "directly" using the out-coupling 
kind of experiment we have been discussing. 



VII. CONCLUSIONS 



This paper has studied the weak coupling of two trapped dilute Bose gases using a tunneling Hamiltonian approach. 
When considering tunneling between two trapped condensates, we obtained a coherent Josephson-type time-dependent 
contribution given in (|3l|), with a frequency Q = oj — Ae — A/i (see also Ref. j[3)). The expression in (|3l|), which 
corresponds to coherent transfer of atoms between two condensates, is first order in the coupling 7. In addition, we 
found that two types of quasiparticle currents emerge from treating a perturbation of the kind shown in to second 
order: a static, time-independent current given by ( f40| ) and ( |4l| ) involving only regular correlation functions of the 

form (ip^ipS, and a time-dependent current given by (^9|) and ( |60[ ) involving a second harmonic with frequency 2f2. 

The latter is associated with anomalous correlation functions of the kind (^t/jip), and is reminiscent of the Josephson 

current carried by Cooper pairs in superconductors [ p~T] | . These tunneling currents were written in terms of single- 
particle spectral density functions that contain all the details of the microscopic correlations of the Bose-condensed 
gases. In this manner, the tunneling dynamics is cleanly separated from the specific approximation used for the 
many-body description of the Bose-condensed gas. This separation was one of the main goals of our work. 

Our paper further serves to unify and extend the work of Rcf. |li| on Josephson-type and other kinds of tunneling 
phenomena in trapped Bose gases, as well as recent work || [|] dealing with out-coupling from a single Bose-condensed 
gas. It provides a transparent way of expressing the various contributions to the tunneling processes involving an 
inhomogeneous trapped Bose gas, before any specific approximation is introduced. The present paper thus provides 
a foundation upon which one can study the details of different tunneling processes within various many-body approx- 
imations for the dynamics of a Bose gas. Our work re-emphasizes the potential usefulness of the kind of out-coupling 
experiments discussed in Refs. m as a probe of the many-body physics of a Bose-condensed gas, including systems 
of reduced dimensionality. 
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As we have noted in the text, our present study contrasts with the work reported in Refs. ||, ||, |9|. The latter 
deal with the density fluctuation spectrum while our present paper (see also Ref. ||, ||) probes the single-particle 
Green's function of a trapped Bose gas, which is perhaps a more basic measure of the condensate-induced correlations 
between atoms. 

For illustration, we have used our formalism in Sections ^ and ^ to find a general expression for the out-coupling 
atom current using the Bogoliubov-Popov quasiparticle approximation. Our results exhibit the same three processes 
obtained by the Oxford group ]3], |I| . The first process is the tunneling of an atom out of the condensate of the trapped 
gas. The second process is the "quantum evaporation" of an atom from the trapped gas, destroying a quasiparticle 
excitation in the trapped gas and creating an atom in the out-coupled gas. The third process, referred to as a "pair- 
breaking" process in Ref. || , involves an atom tunneling out of the trapped gas concurrently with the creation of a 
quasiparticle excitation in the trapped gas (see also Ref. [p5|). 

In Section VI , we have numerically evaluated the relative contributions of these three processes as a function of the 
laser detuning parameter 5 = lu — Ae = cu a — loi, — As, for various values of the wavevector difference between the two 
laser beams q = k a k(, and the temperature T. This kind of experiment would be especially promising as a probe 
of the correlations in a Bose gas in the case of a large scattering length a made possible by working near a Feshbach 
resonance. This enhances the tunneling contribution from the noncondensate atoms and moreover requires the use 
of an improved single-particle correlation function which includes quantum depletion of the condensate (such as the 
second order Beliaev approximation). 

The tunneling currents from the noncondensate atoms are admittedly small and will require high precision exper- 
iments. We argue that the unique information that can be obtained justifies the effort. In future work, we plan on 
using the formalism set up in this paper to discuss the case when the trapped Bose condensate is in a vortex state. 
In addition, information about the single-particle correlation functions in a one dimensional Bose gas pl| is already 
present in the out-coupled tunneling current. 
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FIG. 1: Raman process coupling two atomic hyperfine levels that may or may not be Bose-condensed. 
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FIG. 2: Out-coupling atom current for a ~ 5nm and g = 2/ano- The detuned frequency 5 is defined in (P4[). The plot on 
the left is at T = 0. The plot on the right is at T = 0.9T B ec- In Figs. [||, we take N = 10 6 atoms, lo = 2tt x 60 #2 and 
7 = 0.2wo- 
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FIG. 3: Out-coupling atom current for a ~ 5nm and q = 2/oho- The plot on the left shows each component contributing to 
the out-coupling atom current. The plot on the right shows the total out-coupling atom current. In both plots T = 0.2Tbec- 
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FIG. 4: Out-coupling current for a ~ 5nm and g = 10/a/fo- The plot on the left is at T = 0. The plot on the right is at 
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FIG. 5: Out-coupling atom current for a ~ 5nm and q = W/ciho- The plot on the left shows each component contributing to 
the out-coupling atom current. The plot on the right shows the total out-coupling atom current. In both plots T = 0.2Tbec- 



